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SYSTEM FOR PHASE, GAIN, AND DC OFFSET ERROR CORRECTION FOR A 

QUADRATURE MODULATOR 

FIELD OF THE INVENTION 
[0001] The present invention relates to radio communications technology. More 

particularly, the present invention relates to quadrature modulators and quadrature signal processing 
circuits and systems. 

BACKGROUND OF THE INVENTION 
[0002] Quadrature modulators operate on separate in-phase (I or i) and quadrature (Q or q) 

data corresponding to a single transmit or receive signal. Modern transmitters, such as direct 
conversion radio transmitters and quadrature modulators, are driven by digital circuits at very high 
sampling rates. Even though the in-phase and quadrature data is created digitally in perfect 
quadrature from the same signal, it is distorted by the analog circuits during transmission due to 
phase, gain and direct current (dc or DC) errors inherent in such analog circuits. As a result, 
compensation for these errors must be substantially continuous due to circuit variations over time, 
temperature and ever widening bandwidths. Therefore, an improved method of compensating for 
phase, gain and dc errors associated with the analog circuits of analog transmitters would be a 
significant improvement in the art. 

SUMMARY OF THE INVENTION 

[0003] A quadrature modulator compensation system for compensating the transmission of 

a source signal includes a transmitter which uses a transmitter input signal to generate a transmitter 

output. Calibration circuitry generates a phase error estimate of the transmitter as a function of an 

angle of intersection between an ideal transmitter envelope and an actual transmitter envelope. The 

calibration circuitry also generates a gain error estimate of the transmitter as a function of variation 

in the actual transmitter envelope, as determined by semi-major and semi-minor axes of an elliptical 
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transmitter waveform. Further, the calibration circuitry determines a centroid of the actual 
transmitter envelope, and estimates dc offsets in an in-phase component and a quadrature 
component of the source signal as a function of the centroid of the actual transmitter envelope. 
Predistortion circuitry receives the phase, gain, and dc offset error estimates and the source signal as 
inputs and provides the transmitter input signal as a function of the phase, gain, and dc offset error 
estimates. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0004] FIG. 1 is a block diagram of a quadrature modulator compensation system 

incorporating concepts of the present invention. 

[0005] FIG. 2 is a circuit diagram illustrating one embodiment of the predistortion or 

correction circuitry shown in the system of FIG. 1. 

[0006] FIG. 3 is a circuit diagram illustrating one embodiment of the error transmitter 

shown in the system of FIG. 1 . 

[0007] FIG. 4 is a block diagram illustrating one embodiment of the calibration circuitry 

shown in the system of FIG. 1 . 

[0008] FIG. 5 is a block diagram illustrating one embodiment of circuitry configured to 

implement the find slope function or circuitry shown in FIG. 4. 

[0009] FIG. 6 is a diagrammatic plot illustrating the distortion caused by the phase error 

which makes the transmitter signal elliptical and which can be used to determine phase error 
using the angle of intersection methods of the present invention. 

[0010] FIG. 7 is a block diagram illustrating one embodiment of circuitry configured to 

implement the phase gain computation function or circuitry shown in FIG. 4. 



DETAILED DESCRIPTION OF AN ILLUSTRATIVE EMBODIMENT 
[0011] The present invention includes a method to compute the phase, gain, and dc error 

estimates associated with the analog circuits of quadrature transmitters. The present invention 
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includes, and is described with reference to, an apparatus configured to implement algorithms 
that detect and compensate for these phase, gain, and dc offset error estimates. 
[0012] A general system or block diagram is provided in FIG. 1 in order to illustrate the 

methods of the invention. Complex signals are assumed at the output of each functional block 
shown in FIG. 1, unless otherwise noted in the discussion that follows. System 100 includes 
signal source 110, predistortion or correction circuitry 120, error transmitter 130, ideal receiver 
140, low pass filter 150, envelope detector 160, low pass filter 170, and calibration circuitry 180. 
Each of these functions or functional blocks can be implemented in circuitry or in one or more 
suitably programmed processors, application specific integrated circuits (ASICs), or other 
circuitry. Also, more than one of these functions can be implemented in the same processor or 
other circuitry. References to processors or circuitry being adapted to implement these functions 
means that the processors or circuitry are programmed or configured to implant the functions. 
[0013] Signal source 110 is the data generator that provides data to drive the simulation 

for the following description which explains the methods of the present invention. In the system 
100, signal source 110 provides the signal 111 to be transmitted by transmitter 130 after 
predistortion or correction by circuitry 120. Predistortion circuitry or function 120 changes the 
source data in such a way as to compensate for transmitter errors before transmission occurs. 
Error transmitter 130 adds errors to the signal provided by predistortion or correction circuitry 
120 and translates it to the radio, or channel, frequency (rf) for transmission. Envelope detector 
160 provides feedback about the transmitter's envelope to calibration circuitry or function 180 via 
low pass filter 170 which attenuates double channel and harmonic frequencies. The calibration 
algorithms or methods of the present invention implemented by calibration circuitry 180 compute 
correction coefficients and values which drive a correction model. Ideal receiver 140 and its 
associated low pass filter 150 are included for visual performance results and analysis. Detailed 
explanations of the functions performed by each of these blocks are provided as follows. 
[0014] Signal source 110 illustrated in FIG. 1 is typically an integrated circuit (IC) 

microprocessor optimized for digital signal processing (DSP). Signal source 110 can also be a 
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digital IC such as an ASIC. The complex output s[n] (designated by reference number 1 1 1) of 
source 110 must be predistorted by the second function in FIG. 1 (circuitry 120) in an appropriate 
way before application to error transmitter or transmitter function 130. Signal source 110 also 
provides at output 1 12 its phase Gn to the calibration circuitry or function 180, which is discussed 
later in greater detail For the purposes of discussion, signal s[n] can be described as a complex 
sinusoid, where a> is the modulation frequency, n is the sample number, and Ts is the sampling 
period, using the relationship shown in Equation 1 : 

Equation 1 

s n =cos(oo n-Ts)-hj sin(co n-Ts) 

[0015] The purpose of the predistortion function implemented in block 120 is to correct 

for phase, gain, and dc offset errors inherent in transmitter 130. FIG. 2 is a circuit diagram 
illustrating the predistortion or correction function applied to input signal s[n] by circuitry 120 in 
embodiments of the present invention. As can be seen in FIG. 1, in one embodiment, calibration 
circuitry 180 provides to predistortion circuitry 120 the coefficients or values phEstfk] 181, 
GnEst[k] 182, idcEst[k] 183 and qdcEst[k] 184. Coefficients phEst[k] and GnEst[k] are 
estimates of the phase error and the gain error, respectively, between the i and q quadrature 
components 201 and 202 (FIG. 2) of the input signal s[n] 111. Coefficients idcEst[k] and 
qdcEst[k] are respectively estimates of the dc offset in the i and q components of input signal 
s[n]. The coefficients provide an approximation of the inverse error function of transmitter 130. 
That is, given that the ideal transmitter function output as a function of its input is represented by 
Ti (where Ti = T* input), and given that the error in a practical transmitter is represented by E, 

then the predistortion function 120 approximates E'l in the limit as time goes to infinity. 
Oversimplifying for the moment, from the above assuming that the practical (or error) transmitter 
output is given by Tp (where Tp=T*E*input), and assuming that the predistortion function output 
121 is given by s c [n] (where s c [n] = E"l*s[n]), then the error transmitter output sr^n] 



(designated by reference number 131 in FIG. 1) is given by Equation 2: 

sr d [n] = T*E*E-!*s[n] = T*s[n] 
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Equation 2 



This is the desired result after the appropriate substitutions from FIG. 1 have been made, 
assuming perfect correction. 

[0016] FIG. 2 illustrates the model assumed in predistortion circuitry 120 for the inverse 

error function, E~l, of error transmitter 130 described above. This model results from assuming 
that the phase and gain errors are modeled in the quadrature leg of the transmitter as shown in 
Figure 3. The inputs i,q (201,202) and outputs i c? q c (203,204) represent quadrature components 
of input signal s[n] and output signal s c [n], respectively, in FIG. 1. Kl is a gain error term given 

by Equation 3: 

Equation 3 

Kl^GnEst^cos fphEsO 



K2 is a crossover term involved in phase error given by Equation 4: 



K2 k =GnEst k sin f phEst R 



Equation 4 



S is a gain correction term to maintain unity gain given by Equation 5: 



S k = f GnEsycos f phEst k 



-l 



Equation 5 



Assuming that s =i c +j *q .and assuming the above expressions for coefficients Kl, K2, and 

n n n 

S in FIG. 2, the output 203,204 (S c [n] 121 in FIG. 1) is determined using the relationship of 
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Equation 6: 

Equation 6 

s c =S (Kl i n + K2-q n ) - idcEst^j -(s.q n - qdcEst k ) 



[0017] The estimates of the offset errors idcEst and qdcEst are independent of the phase 

and gain errors and are computed as the centroid of the received envelope (from envelope detector 
160). As will be described below in greater detail, the envelope of error transmitter 130 is resolved 
to the signal source's phase with delay compensation in order to approximate the errors as shown 
here* 

[0018] Error transmitter 130 implements an error transmitter function which can be 

modeled mathematically as shown in FIG. 3. This function translates the baseband signal sr<j[n] 
(i.e., component signals i c and q c ) at output 131 to a local oscillator (LO), or channel, frequency 
Go. The transmitter or transmitter function may be complex, but the real model is used in 
forming the real output 131, sr^n], shown in FIG. 1. This model assumes all phase & gain 
errors are in the quadrature channel. That is, the in-phase channel i is used as the reference 
channel All of the models discussed here are inter-related using this assumption. However, in 
view of this disclosure, new models can be derived in order to model the phase and gain errors in 
the in-phase channel. The output srd[n] 131 of error transmitter 130 illustrated in FIGS. 1 and 3 
can be expressed using the relationship shown in Equation 7: 

Equation 7 

sr d [*\ = Qc, +^)cos(<9 0n )-G( 9cw + q dc )$in(0 On +</>) 

where idc and qdc represent the actual in-phase and quadrature dc offset errors, G represents the 
actual gain error, and <[> represents the actual phase error inherent in the transmitter (these are the 
errors to be estimated by calibration circuitry 180). 

[0019] Envelope detector 160 in FIG. 1 is applied to the error transmitter output, srjtn], 
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to provide feedback about the transmitter output 131 to a calibration function or circuitry 180. 
The value G d provided to calibration circuitry 180 is a gain term to compensate for using a real 
signal as opposed to a complex signal. The gain & dc errors, G r and DC r , respectively, inherent 
in envelope detector 160 are factored into its real output, Er[n] 161, which is filtered by low pass 
filter 170 to remove any double frequency and harmonic terms before use by the Calibration 
function as Efr [n] 171. 

[0020] The calibration circuitry function 180 shown in FIG. 1 is also shown separately in 

FIG. 4. This circuitry resolves the filtered envelope detector output, Efr[n] 171, into a quadrature 
signal, compensates for circuit and filter delays, and computes updated estimates of the phase 
(phEst[k]), gain (GnEst[k]), and dc offset (idcEst[k] & qdcEst[k]) errors at a decimated rate. The 
phase error detection and correction algorithm implemented by circuitry 180 computes the angle 
of intersection of the envelope of the transmitter output and the ideal transmitter envelope which 
is then integrated to form a phase error estimate. The gain error detection and correction 
algorithm implemented by circuitry 180 computes the maximum envelope variation of the 
transmitter output which is then integrated to form a gain error estimate. The dc offset error 
detection and correction algorithm implemented by circuitry 180 computes the centroid of the 
envelope of the transmitter output which is then integrated to form a dc offset error estimate for 
both the i and q channels. The calibration circuitry or function 180 implements several simple 
functions and a discussion of these individual components is provided below with reference to 
FIG. 4. FIG. 4 shows a functional decomposition of circuitry 180 and is the top level diagram for 
purposes of discussion. 

[0021] FIG. 4 illustrates various functions implemented by calibration circuitry 180 in 

accordance with aspects of the present invention. The various functional blocks illustrated in 
FIG. 4 can be implemented within one or more processors in accordance with the following 
discussion. However, in other embodiments, each functional block is implemented using 
separate circuitry. While all such embodiments are considered to be within the scope of the 
present invention, the following discussion refers to the separate functions illustrated in FIG. 4 as 
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being separate circuitry. 

[0022] As can be seen in FIG. 4, calibration circuitry 180 includes polar to rectangular 

converter 405. This circuitry or function uses the phase 9[n] from signal source 110 and a 
constant phase delay 9dly to resolve the filtered envelope Efr[n] into quadrature components 
needed by the other calibration function elements. The phase delay Odly represents hardware 
delay and can be retrieved from a table of values stored in memory which are accessed by circuit 
405. The result, s r [n] provided by polar to rectangular converter 405 at output 406 includes 

projection of the filtered envelope phasor onto the real and imaginary axes. A mathematical 
model for this process is shown in Equation 8: 

Equation 8 

s r =Efycos ^9 n + 0dly j + j -Efi-sin (e n + Odly) 

[0023] Sample delay circuitry 410 receives signal s[n] from signal source 1 10, as well as 

delay, dly, to compensate for pure time (i.e., sample) delay when needed. Dly is a pure time delay 
representing the time compensation required to allow the association of the proper samples when 
computing coordinates of intersection. It is an artifact of hardware and filtering delays from the 
signal source to the calibration function that can be determined off-line and stored in a table or 
constant. The output ss[n] 411 can be defined for example using the relationship shown in 
Equation 9. 

Equation 9 

ss =s „ 

n n— dly 

[0024] Envelope detector 415 is used to detect the quadrature signal S r [n] 406 formed 

during the resolution of the filtered envelope Efr[n] 171 to the source's phase and compensation 
for delay. This function may not be needed depending on sampling resolution, timing and phase 
delays. The real and imaginary operators are a modeling convenience. Any implementation will 
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normally use rectangular coordinates. The envelope output 416, Env r [n], of envelope detector 
415 is given by: 

Equation 10 

where, 

Re(Sm) is the real portion of S r [n]; and 
Im(Sm) is the imaginary portion of S r [n]. 

[0025] Envelope detector 420 is used to compute the envelope of source signal S[n] 111, 

including any necessary sample delay 0dly, which will be used as the desired or reference 
waveform. The mathematical model can be the same as shown in Equation 10, with the 
appropriate substitutions made at the input and output (ss[n] for S r [n], and Env t [n] for Env r [n] 
respectively). 

[0026] Intersection identifying circuitry 425 finds the points of intersection of the 

received and filtered transmitter envelope, Env r [n] 416, with the desired waveform, Envtfn] 421. 

Index of intersections ix (ix[k] in FIG. 4) can be calculated as shown in Equation 1 1 below. The 
u() function in Equation 11 is the unit step function used to detect when the envelope error 
variable, Env e [n], changes sign. The envelope error variable Env e [n] is defined in Equation 12. 

The input to u() is always positive except when there is a sample on each side of an intersection, 
in which case it is negative. Whenever the error changes sign, an intersection is indicated, and 
the index of that intersection, n, is saved in ix 426. 

Equation 1 1 

ix = -n - [u{Env en - Env e ^ )- 1) 
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Equation 12 

Env^ = Env „ _ Env+. 



[0027] Slope identifying circuitry 430 finds the slope of each intersecting curve s r [n] and 

ss[n] at the point of intersection (n=ix). Any intersecting point in one phasor revolution is 
sufficient. This method uses the intersecting point sample S r [ix] (intersecting point in curve 
S r [n]) and SS[ix-l] (intersecting point in curve SS[n]) and the previous sample in the 
approximation of the slope, but linear interpolation could be used to find two points on either 
side of the intersecting point on the respective curves. This latter method could be used to give a 
more accurate slope value at the point of intersection between the curves and to possibly reduce 
the required sampling rate. The complex outputs mrc[ix] 432 and mtc[ix] 433 of slope 
identifying function 430 can be determined using the relationships of Equations 13 and 14, which 
are modeled in FIG. 5. 

Equation 13 

mrc=Re/s r ^-Refs r ) + \ -(lm(s r \-Ws r )) 
" I r ixj { r ix-ij [ [ r ixj \ r ix-i I) 

where, 

Sn x =Sr[n] for n=ix (point of intersection) and; 
Srix-i=S r [n] for n=ix-l (previous sample). 



Equation 14 



mtc. =Re ss. -Re ss. J+j • Imfss. -Im ss. 

IX \ ixj V ix— 1/ J \ \ IX/ \ ix— 1 



where, 

SSix^SSfn] for n=ix; and 
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SSix-i=SS[n] for n=ix-l (previous sample). 

[0028] The phase error of the transmitter expresses itself in the angle of intersection of 

the desired and actual transmitter envelopes Envt[n] and Env r [n], respectively. As shown in FIG. 
6, the distortion caused by the phase error makes the transmitter signal elliptical with a 45 degree 
tilt (see ellipse 605). The angle of intersection, a between the transmitter signal (represented by 
605) and the desired transmitter signal (represented by 610) can be computed from any point of 
intersection of the two curves as derived from plane geometry discussed below. 
The sum of all the angles of triangle 615 formed by angles a, pi and p3 is 180 degrees. Angle 
p3 is equal to 180-p2, and the angle of intersection, a, is therefore equal to p2-pl . 
[0029] It can be shown that: 

Equation 15 

tan( a )=tan( p 2 - p 1 )=( tan( p 2) - tan( p 1 ) ) •( 1 + tan( p 2) tan( p 1 ))"* 



Now, the tangent value effectively gives the slope of the line 620 through the point of 
intersection 625, and we will let mr represent tan(p2) & mt represent tan(pl) giving: 

Equation 16 



tan 



- 1 / mr- mt \ 
—__ i 

l + mrmt/ 



•INV 



[0030] Where INV is discussed below and where the slopes are formed as shown in 

Equations 17 and 18 (division by zero must be avoided). 

Equation 17 



ml k= Im ( mtC ix)-( Re ( mtC ix)) 
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Equation 18 



mi k =Im(mrcx x )fRe(mrc ix ) 1 " 1 



[0031] Therefore, assuming that the angle of intersection a of the two waveforms 

approximates the transmitter phase error, then the slope of the intersections can be computed 
using the above equation thus giving a method of estimating the phase error. The INV factor in 
the Equation 16 is needed because there is a phase ambiguity depending upon which intersection 
point is used. If the magnitude of the imaginary coordinate of the point of intersection 625 is 
greater than the magnitude of the real part, INV = -1, otherwise it is 1 . 

[0032] A circuit 435 for implementing the mathematics described above is shown in 

FIGS. 4 and 7. Since the signum function outputs a -1 for a negative input, 0 for a 0 input, and 1 
for a positive input, the impulse function is included to make 0 inputs positive. Each update of 
ad is integrated with the previous phase error estimate, thus forming a new one. The phase 
estimate phEst[k] 181 computed in this way converges on the actual phase error to within 0.1 
degrees. The phase estimate phEst[k] so computed will be filtered (not shown) to provide noise 
immunity. 

[0033] Referring back to FIG. 4, find max and min function or circuitry 440 is included 

in circuitry 180, receiving as an input transmitter envelope Env r [n], Circuitry 440 finds the 
variation in the transmitter envelope for use in approximating the gain error in the error 
transmitter 130. It compares the filtered envelope with a maximum and a minimum value to 
effectively compute the semi-major and semi-minor axes (Ex[k] 441 and En[k] 442 respectively) 
of the elliptical transmitter waveform. In computing Ex, the unit step function is 1 for a 0 or 
positive input (that is, when Env r [n] > Ex[k]) which selects the nth value of the envelope to 
replace Ex[k]. There is no change in Ex[k] when Env r [n] < Ex[k]. 

[0034] Calibration circuitry 180 also includes gain estimate computation circuitry 445. 
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Circuitry 445 receives as inputs from circuitry 440 the semi-major and semi-minor axes Ex[k] 
441 and En[k] 442. Circuit 445 integrates the envelope variation with the old gain estimate after 
the phase is toggled for each iteration. The phase inversion which may be replaced by a 1-pole 
averager is represented as EvenOdd which alternately inverts the sign of the envelope variation 
and is needed to maintain stability in the gain estimates because the phase estimate can be 
occuring simultaneously and the two estimates interact. The gain estimate GnEstfk] 182 provided 
by circuitry 445 can converge to within 0.01 db of the actual gain G. It can be filtered (not 
shown) to provide noise immunity similar to the handling of the Phase Estimate mentioned 
above. The output is given by Equation 19: 

Equation 19 

GnEst k =GnEst k _ j + EvenOdd • (E^ - EnJ 

[0035] Also included in calibration circuitry 180 are centroid computation circuitry 450 

and dc estimate generating circuitry 455. Circuitry 450 implements a mathematical model for 
computing the centroid of the filtered transmitter envelope. The components of the Calculus of 
the centroid of a curve include the length, L, of the curve and the moments of each element of L 
with respect to each of two aces (I and Q in this case). The moment of each element of the 
waveform is computed for both the in-phase and quadrature components of the envelope. The 
length L, of the path traced out by the received and filtered, transmitter envelope phasor, s,- (406), 
is also computed. L is used in the next two computations of the centroid of the phasor path, idc 
and qdc. Other methods of effectively computing the centroid by using the envelope error, Env e , 
may be possible. The resulting updates are then integrated by the dc estimate computation 
circuitry 455 with the previous dc offset estimates for use during the next epoch or decimation 
period. The envelope notation used here, s r (406), results from the quadrature resolution of 
Efr[n] 171 . The number of samples in a rotation of the envelope phasor is 0.5 *n x where n x is the 
total number of samples generated for study (two rotations, or 4*pi, worth of samples). The last 
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half (thus the 0,5 factor on n x below) are used to avoid transient effects in the data. 
[0036] The waveform path length L is then given by Equation 20: 

Equation 20 

n x - 1 

v 2 J( ta ( s 'J-""( s ' b -,)) 2+ ( Re ( s vi- Re ( s ' k -)) 2 

b-0.5n x 

where nx is defined above and s r [b] and s r [b-l] (406) are adjacent samples of the quadrature 
resolution of Efr (171). 

[0037] The quadrature components of the centroid of the received and filtered envelope 

waveform can then be computed as shown in Equations 21 and 22 which effectively satisfy the 
evaluation of the well known, however more recognizable, equations for the centrod of a curve in 
an x-y plane (x c = Jx*ds and y c = Jy*ds) in discrete form : 



Equation 21 



n x -i 



b = 0.5-n , 



{ V \ T b-i}j 



Ws 



n 



r by \ r b-i 



+ (Re[s 



-Re/s r 



\2 



■w 



v 1 



v- 2 



b = 0.5n. 



Ws 

I 



V +fa \'Vi. 



•0.5 [Ws r \- Ws r ^ 2 



> Rqs 



-Re's 



\ b-i 



Equation 22 



[0038] The centroid coordinates are integrated with the appropriate previous dc error 

estimates. The dc estimates also will be filtered (not shown) to provide noise immunity similar 
to the handling of the phase and gain estimates mentioned above. The outputs 183 and 184 can 
be given by Equations 23 and 24; 
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Equation 23 

idc est = idc est t + idc k 

k k— 1 

Equation 24 

^dc estk =qdc estki +qdc k 

[0039] The invention has been described as generating a phase error estimate, a gain error 

estimate, and a dc offset estimate. Said estimates may be generated simultaneously or in any order. 
Furthermore, any combination of said estimates may be performed as desired. 
[0040] Although the present invention has been described with reference to illustrative 

embodiments, workers skilled in the art will recognize that changes may be made in form and detail 
without departing from the spirit and scope of the invention. 
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